A Change of Rotation Profile in the Envelope in the HH 111 
Protostellar System: A Transition to a Disk? 
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ABSTRACT 

The HH 111 protostellar system consists of two Class I sources (VLA 1 and 2) 
with putative disks deeply embedded in a flattened envelope at a distance of 400 
pc. Here is a follow-up study of this system in C 18 (J=2-l), SO (Nj = 5 6 — 4 5 ), 
and 1.33 mm continuum at ~ 1" (400 AU) resolution, and it may show for the 
first time how a rotationally supported disk can be formed inside an infalling 
envelope. The 1.33 mm continuum emission is seen arisen from both sources, 
likely tracing the dusty putative disks around them. In particular, the emission 
around the VLA 1 source is elongated in the equatorial plane with a radius of ~ 
300 AU. The envelope is well seen in C 18 0, extending to ~ 7000 AU out from the 
VLA 1 source, with the innermost part overlapping with the dusty disk. It has 
a differential rotation, with the outer part (~ 2000-7000 AU) better described 
by a rotation that has constant specific angular momentum and the inner part 
(~ 60-2000 AU) by a Keplerian rotation. The envelope seems to also have some 
infall motion that is smaller than the rotation motion. Thus, the material in the 
outer part of the envelope seems to be slowly spiraling inward with its angular 
momentum and the rotation can indeed become Keplerian in the inner part. A 
compact SO emission is seen around the VLA 1 source with a radius of ~ 400 
AU and it may trace a shock such as an (inner) accretion shock around the disk. 

Subject headings: circumstellar matter - stars: formation — ISM: individual (HH 
111) 



1. Introduction 



Stars are formed inside molecular cloud cores by means of gravitational collapse. The 
details of the process, however, are complicated by the presence of magnetic fields and angular 
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momentum. As a result, in addition to infall (or collapse), rotation and outflow are also seen 
toward star-forming regions. In theory, a rotationally supported disk is expected to form 
inside a collapsing core around a protostar, from which part of the material is accreted by the 
protostar and part is ejected away. Observationally, however, when and how a rotationally 
supported disk is actually formed are still unclear, because of the lack of detailed kinematic 
studies inside the collapsing core. 

This paper is a follow-up stud y of the HH 111 protostellar system. The properties of 
this system have been reported in iPaper II and only the important ones are summarized 
here. This system is deeply embedded in a compact molecular cloud core in the L1617 
cloud of Orion at a distance of 400 pc. It is a Class I protostellar system with a flattened 
envelope, putative disks, and highly collimated jets. At the center of this system, there are 



two sources, VLA 1 and VL A 2, with a pro 



former driving the HH 111 jet (IReipurth et al- 



erted separation of ~ 3" (1200 AU) and the 



19991 ) . Previous observation of this system has 



shown that the flattened envelope seems to be in transition to a rotationally supported disk 
near the VLA 1 source (Lee et al 2009, hereafter Paper I). Here is a follow-up study of this 
system with more complete uv coverage, extending to both shorter and longer uv spacings, 
in 1.33 mm continuum, C 18 (J = 2 — 1), and S O (Nj = 5$ — 4 5 ) emission obtained with 
the Submillimeter Array (SMA)0 ( Ho et al.l 120041 ) . This study may show for the first time 
a change of rotation profile in an envelope and thus how a disk can be formed inside an 
envelope. 



Observations 



Observation toward the HH 111 system was first carried out with t he SMA on 2005 
December 5 in the compact-north configuration, with the results reported in lPaper 1 Follow- 
up observations were then carried out on 2010 January 24 in the extended configuration for 
longer uv spacings to resolve the structure of the envelope and on 2010 March 26 in the 
sub-compact configuration for shorter uv spacings to restore the large-scale structure. The 
230 GHz band receivers were used to observe the 12 CO ( J = 2 - 1), 13 CO ( J = 2 - 1), C 18 
( J = 2 — 1), and SO (Nj = 5q — 4 5 ) lines simultaneously with t he 1.33 mm c ontinuum. The 
rest frequency of the SO line is assumed to be 219.949391 GHz ( )Lovaslll992l ). Note that the 
12 CO and 13 CO lines trace mainly the outflow interaction and will be presented in a future 
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publication. 

In the follow-up observations, the velocity resolution in SO line was increased to ~ 0.28 
km s -1 per channel, similar to that of C 18 0. Also, only one single pointing toward the central 
region was observed in order to have enough sensitivity. The visibility data were calibrated 
with the MIR package. The flux uncertainty was estimated to be ~ 15%. The calibrated 
visibility data were combined with that in t he compa ct-north configuration and then imaged 



with the MIRIAD package, as described in iPaper J . With various different weightings, the 



synthesized beams can have a size from l!'6xl"5 to l''2xl''0. The rms noise levels are ~ 
0.06 Jy beam" 1 in the C 18 and SO channel maps with a beam of l''3xl''l, and 1.6 mJy 
beam -1 in the continuum map with a beam of l''4xl''l. The velocities in the channel maps 
are LSR. 



Results 



As in IPaper II . the results are presented in comparison to a mosaic image based on 
the Hubble Space Telesc ope (HST) NICMOS image ([Fell] 1.64 fim + H 2 at 2.12 /im + 
continuum) obtained by iReipurth et al.l (119991 ). which shows two sources in the infrared, 
reflection nebulae that trace the illuminated outflow cavity walls, and the jet in the system. 
The two sources were also detected in 3.6 cm continuum by the VLA as the VLA 1 and 
2 sources, respectively, at a( 2 ooo) = 05 h 51 m 46f254, foannn) = +02°48'29"65 and «( 2 ooo) = 
05 h 51 m 46?07, (5 (2 ooo) = +02°48'30"76 faeipurth et al.l ll999l : [Rodriguez et all l2008h . These 



VLA positions are more accurate than the NICMOS positions and are thus adopted here as 
the source positions. The systemic velocity in this system is assumed to be 8.90±0.14 km 
s -1 LSR, as discussed in Section [3.2. II Throughout this paper, the velocity is relative to this 
systemic value. 



3.1. 1.33 mm Continuum Emission 

As shown in Figure [TJ the continuum emission is now spatially resolved into two compo- 
nents, a bright one associated with the VLA 1 source and a much fainter one (~5a detection 
with 1 cr = 1.6 mJy beam -1 ) associated with the VLA 2 source. The one associated with the 
VLA 1 source is slightly elongated perpendicular to the jet axis. This elongation is better 
seen in the CLEAN component map, which shows a faint disk-like structure extending to 
~ 0"7 (280 AU) to the north and 0"9 (360 AU) to the south. Note that in the CLEAN 
component map, the central peak at the VLA 1 position has a FWHM size of < 0!'6 (240 
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AU) and ~ 90% of the flux, but no emission is detected toward the dark ridge position at 
the cent er of the reflection nebulae likely because of a lack of sensitivity. As discussed in 



Paper t . the emission is thermal emission from dust. Thus, the emission is believed to come 
mainly from a dus ty disk, with the inne r part of which already seen with a size of ~ 0''15 
(60 AU) at 7 mm ([Rodriguez et all 120081) . The total flux of this emission is ~ 300±40 mJy, 
sli ghtly hig her than that found in lPaper I] because of the inclusion of the shorter uv spacings. 
In iPaper t . the dust was found to have a temperature of 41-64 K, with a mass opacity of 
0.026 cm 2 g -1 and the emission being optically thin at 1.33 mm. Therefore, the gas and 
dust associated with this emission have a mass of ~ 0.13 ± 0.03 M & . Most of this mass, 
however, is from the central unresolved peak. On the other hand, the total flux toward the 
VLA 2 source is ~ 11±2 mJy. The mass would be ~ 0.005 M Q if the emission has the same 
properties as that around the VLA 1 source. 



3.2. Line Emission 

3.2.1. Spectra 

Figure [2] shows the C 18 and SO spectra toward the VLA 1 source. These lines are 
detected upto ~ ±4.0 km s _1 from the systemic velocity. The SO spectrum is rather sym- 
metric and can be roughly described by a Gaussian line profile, and its line center can be 
used here to define the systemic velocity. Since part of the SO emission comes from the 
outflow as discussed later, further observation at higher resolution is needed to refine the 
systemic velocity by separating the outflow from the envelope. The C 18 spectrum shows a 
double-peaked profile with a dip near the systemic velocity, and with the blu eshifted p eak 



brighter than the redshifted peak (i.e., the blue asymmetry). As discussed in IPaper II . the 
dip with the blue asymmetry could be due to self-absorption because of an infall motion in 
the envelope and the missing flux resolved out by the interferometer. The SO emission does 
not show a dip, likely because it is optically thin due to low SO abundance. 



3.2.2. Morphologies at different velocities 

In order to see how the envelope structure changes with velocity, the line emission is 
divided into 3 velocity ranges: low (0—1 km s _1 ), medium (1—2 km s _1 ), and high (2—3.7 
km s -1 ) velocity ranges, on the redshifted and blueshifted sides. 

The structure of the envelope can be clearly seen in C 18 at low velocity, as shown 
in Figure |3k. The envelope is seen extending to ~ 16" (6400 AU) to the north and south 
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from the VLA 1 source, about twice as extended as that seen in lPaper t . It is well confined 
in the equatorial plane perpendicular to the jet axis, fitting right in between the cavity 
walls outlined by the reflection nebulae. The detailed structure, however, is complicated. 
In particular, the emission in the north beyond 12" out from the source is deviated away 
from the equatorial plane and the emission in the south is curvy. With respect to the VLA 
1 source, the blueshifted emission is to the north and the redshifted one is to the south, as 
expected for a rotating envelope. As we zoom into the inner region at higher resolution, 
the emission appears to be aligned with the dark ridge at ~ 0"5 in the west of the VLA 
1 source (see Figure [3b and c for the central region). The blueshifted emission also arises 
from the outflow, extending to the west along the jet axis, to the northeast along the cavity 
wall, and to the southwest. At medium velocity, the emission shrinks toward the VLA 1 
source (Figure EH), as expected if the velocity of the emission increases toward the source. 
The blueshifted peak is closer to the source than the redshifted peak. In addition, both the 
blueshifted and redshifted peaks shift to the east toward the equatorial plane of the VLA 1 
source, as compared to those at low velocity. Note that the blueshifted emission also extends 
toward the VLA 2 source and along the jet. At high velocity, the emission shrinks even closer 
to the source (Figure [3k). In addition, the blueshifted and redshifted peaks shift further to 
the east into the equatorial plane of the VLA 1 source, coincident with the disk-like structure 
seen in the continuum in the CLEAN component map. Thus, the motion at high velocity 
is highly dominated by the rotation, with the emission in the equatorial plane. Again, the 
blueshifted peak is closer to the source than the redshifted peak, and here by ~ 0''3. 

Unlike the C 18 emission, the SO emission is compact with a size of < 3", tracing the 
innermost part of the envelope. At low velocity, the blueshifted and redshifted peaks are 
seen on the opposite sides of the VLA 1 source in the equatorial plane (Figure [3f), tracing 
the rotation motion. That both the redshifted and blueshifted emission also extend to the 
opposite sides of the source could be due to an additional motion, either infall or outflow. In 
addition, the blueshifted emission also extends to the VLA 2 source. At medium velocity, the 
emission shrinks toward the VLA 1 source (Figure [3g). The peak of the blueshifted emission, 
however, is shifted to the west of the VLA 1 source, likely due to the contamination of the 
emission extending to the VLA 2 source and the emission at the outflow base. At high 
velocity, the emission is very weak (Figure Eh). Both the redshifted and blueshifted emission 
seem to be contaminated by the outflow emission, extending along the cavity walls and along 
the jet axis. 
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3.2.3. Rotation motion 



The rotation profile in the envelope can now be better constrained than that in lPaper J . 
with more complete uv coverage. Figures UK through Hp show the position-velocity (PV) 
diagrams centered at the VLA 1 source cut perpendicular to the jet axis in C 18 0, with 
increasing resolution toward the inner region. The blueshifted emission is seen mainly to 
the north and redshifted emission mainly to the south, confirming that the motion in the 
envelope is dominated by rotation. The blueshifted emission beyond ~ 12" in the north shows 
a velocity increasing with the distance from the source (Figure HK). That part of the emission 
is seen deviated from the equatorial plane (see Fig. |3K) and thus may trace an interaction 
with the surrounding. Except for that part, the rotation velocity increases toward the source 
and thus could be Keplerian or that has constant s pecific angular momentum as seen in some 



infalfing envelop es in the Class I phase (see, e.g., lOhashi et al.lll997al ; iMomose et al.lll998 



Lee fc Hcll2005n . 



Figures [5K and [5b show the rotation velocity (v^) and the implied specific angular 
momentum (I = Rv^), respectively, as a function of R, which is the radius measured from 
the source. The data points are extracted from the PV diagrams with increasing resolution 
toward the inner region and marked as "x" s in Figures HK through Hp. Both the rotation 
velocity and specific angular momentum are corrected for the inclination angle assuming 
that the envelope is perpendicular to the jet axis and thus has an inclination angle of 10° to 



the line of sight (IReipurth. Raga. fc Heathcotdll992l ). As can be seen, the rotation profile 



can be roughly divided into two parts. For R > 5" (2000 AU), the rotation velocity can be 
roughly reproduced by the rotation that has constant specific angular momentum (see also 
Figures HK through Hp), with ~ 3.90( J R/i? )~ 1 km s" 1 and #0 = 1" (or 400 AU). The 
specific angular momentum there is almost constant at ~ 1550 AU km s _1 . For R < 5", the 
rotation velocity increases slower toward the source and is better described by a Keplerian 
law (see also FiguresHK throughHb), with ~ 1.75(-R/-Ro)~ ' 5 km s _1 . The specific angular 
momentum decreases steadily to ~ 300—400 AU km s _1 at ~ 100 AU (0''25). Note that in 
Figures [5K and [5b, for R < 1", the blueshifted data points and the redshifted data points 
are separated from each other. As mentioned, at high velocity, the peak of the blueshifted 
emission in the north is closer to the source by ~ 0!'3 than that of the redshifted emission in 
the south (Figure [3k). This difference in distance is also seen in the PV diagram in Figure 
Hp. Thus, it is possible that the center of the envelope is actually located at 0''15 to the 
south of the VLA 1 source. If so, then the blueshifted data points and the redshifted data 
points will move closer to each other, as seen in Figures [5b and [5bl. Observations at high 
angular resolution are needed to confirm this. 



Rotation is also seen in the compact SO emission near the source. As shown in Figure 
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HH, the PV structure of the SO emission overlaps with that of the C 18 emission near the 
source, tracing the rotation there. As in C 18 0, the center of the PV structure seems to be 
located slightly to the south. 



3.2.4- Infall motion 

The infall motion can be studied with the PV diagrams cut along the jet axis in C 18 
and SO, as shown in Figure Hfe. The C 18 emission at ~ — 1 km s _1 on the blueshifted 
side with the position offset < —1" is from the outflow emission extending to the west along 
the jet axis as seen in the low and medium- velocity maps (see Figures and d). This 
blueshifted emission in the west also extends close to the source, contaminating the envelope 
emission. Also, the SO emission on the blueshifted side at ~ (—0.6 km s" 1 , — 1''5) extends 
to the VLA 2 source, and thus should be excluded for studying the infall motion. Excluding 
those parts of the emission, the C 18 and SO PV structures appear to be similar and both 
are roughly symmetric about the source, with the blueshifted emission shifted slightly to the 
west and the redshifted to the east. With the near end of the envelope tilted to the east, 
these PV structures, together with the double-peaked line profile with a blue asymmetry and 
an absorption dip at th e systemic velocity (Fig [2]), suggest an infall motion in the envelope, 



as discussed in iPaper I . 



4. A flattened rotating envelope model 

In order to roughly derive the structure, the physical properties (i.e., density and tem- 
perature), and the infall velocity of the envelope, a simple flattened envelope model is used to 
reproduce the PV diagrams, spectrum, and integrated intensity map of the C 18 emission. 
In this model, the envelope has an inner radius of -R in , an outer radius of .Rout, and a half- 
opening angle of 8q, so that its half thickness, H, can be given in the cylindrical coordinates 
(R,9,z) by 

H(R) = max(# Q ,.Rtan0o) (1) 

where 8q is measured from the equatorial plane, and H is the minimum value of the half 
thickness near the source. The number density of molecular hydrogen in the envelope is 
assumed to be given by 

■VW+ 



n(R,z)=n {^-^ J (2) 

where Rq = 400 AU (1", as in Section I3.2.3H . n is the density at Rq and p is a power 



law index assumed to be —1.5, as in many theoretical infalling models (see, e.g., IShulll977 
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Nakamural 120001 ) . The abundance of C O rel ative to molecular hydrogen is assumed to be 



constant and given by 1.7 x 10 7 , as in lPaper J . The temperature of the envelope is uncertain 
and assumed to be given by 

t(R) = U^y (3) 

where T is the temperature at Rq and q is a power-law index. 

As discussed in Section 13. 2. 3[ the rotation is assumed to change from that with constant 
specific angular momentum in the outer part to that of Keplerian in the inner part, i.e., 

f ^(if r ' 5 if r < Rt, 

< 4 » 



with the transition radius 



Rt = Ro(-) 2 (5) 

Vk 

at which the Keplerian rotation is the same as that with constant specific angular momentum. 
The infall motion is uncertain. It is assumed to be in the radial direction and given by 



v r {R, z) = v r0 (6) 



Ri 

as in a collapsing envelope. 

In the model, radiative transfer is used to calculate the emission, with an assumption 
of local thermal equilibrium. For simplicity, the line width is assumed to be given by the 
thermal line width only. The channel maps of the emission derived from the model are used 
to calculate the model visibility data with the observed uv coverage. This model visibility 
data are used to obtain the channel maps, and then the integrated intensity map, spectrum, 
and PV diagrams to be compared with the observations. 

Figure [7] shows the comparison of the best model with the observations. As mentioned 
earlier, the envelope is assumed to have an inclination angle of 10° to the line of sight, with 
the nearside tilted to the east. The parameters that are used to reproduce the observations 
are V/~ ~ 1.75 ± 0.2 km s -1 , v c ~ 3.90 ± 0.4 km s _1 , and t> r o ~ —0.7 ± 0.2 km s -1 for the 
envelope kinematics, R in ~ 0'.'15 (60 AU), R out ~ 18" (7200 AU), H ~ 0?4 (160 AU), and 
9 ~ 8±2° for the envelope structure, and n ~ 1.3 ± 0.3 x 10 7 cm -3 , T ~ 40 K, and 
q ~ —0.65 ±0.1 for the envelope physical properties (see also Figure [6] for the parameter 
profiles). The values of Vk and v c are set to those found in section I3.2.3[ resulting in the 
transition radius of Rt ~ 5" (2000 AU). The mass that supports the Keplerian rotation is 
~ 1.38 M Q . Thus, the mass of the central source, which can be estimated by subtracting 
from this mass the mass of the dusty disk close to the source, is ~ 1.25 M . Here, the 
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inner radius R m is set to where it has a Keplerian rotation velocity of ~ 4.5 km s" 1 . To 
is set to be the lower end of the dust temperature (IPaper J ). As can be seen, the simple 
model can reproduce the observed PV diagrams, spectrum, and integrated intensity map 
reasonably well. In the model, the intensity near the systemic velocity is higher than that 
in the observations, probably because a large-scale envelope (with a size larger than 20") 
or core is needed to resolve out the emission there. Also, the blueshifted emission around 
— 1 km s _1 is less than that in the observations because the observed blueshifted emission 
there is contaminated by the outflow emission (Figure [7b). Moreover, as compared to the 
observations, the PV structure perpendicular to the jet axis is shifted slightly to the north 
(Figure [7b), probably because the center of the envelope is indeed located slightly to the 
south, as discussed earlier. 

As a result, the C 18 envelope can be considered as a flattened envelope that extends 
from ~ 60 to 7000 AU out from the source, with a small opening angle of ~ 16°. It has a 
differential rotation, with the outer part (~ 2000-7000 AU) better described by a rotation 
that has constant specific angular momentum and the inner part (~ 60-2000 AU) by a 
Keplerian rotation. The transition radius between the two parts is ~ 2000 AU. The mass of 
the envelope (including helium and molecular hydrogen) is given by the following integration: 

f R f H \/B 2 + z 2 
M e = 2.8m H n / / ( V )~ L5 2nRdRdz 

JR m J-H Rq 

Thus, the mass of the whole envelope extending from 60 to 7000 AU is ~ 0.62 M , about 
45% of the mass of the central source. On the other hand, the mass of the inner part of the 
envelope extending from 60 to 2000 AU is ~ 0.13 M Q . Since this mass is only ~ 10% of the 
mass of the central source, the rotation in the inner part of the envelope can indeed become 
Keplerian. The envelope seems to have a small infall motion. The amount of the infall 
motion, however, could be overestimated due to the outflow contamination in the unresolved 
PV structure along the jet axis. Observation at higher resolution is really needed to confirm 
the infall motion by separating the outflow from the envelope. 



Discussion 



5.1. Comparing to the previous results 



With a more complete uv coverage, the envelope is detecte d to a two times larger 
extent at upto three times higher resolution than that in IPaper t . The infall velocity can 
now be better determined with less outflow contamination, and is found to be ~ 40% of that 



derived in IPaper J . The rotation profile in the inner part of the envelope can now be better 
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constrained with a larger extent of the envelo pe, and i s found to be more like Keplerian but 
with a 30% higher velocity than that found in iPaper t . Thus, the inner part of the envelope 
is no longer to b e consist ent with a dynamically infalling disklike envelope (or pseudo-disk), 
as suggested in IPaper J . Previously, the mass of the central source was derived from the 
infall velocity assuming a dynamical infall and found to be ~ 0.8 M Q , but now is from 
the Keplerian rotation velocity and found to be ~ 1.25 M . In addition, the compact SO 
envelope near the source is now better resolved and thus better separated from the outflow 
cavity walls. The PV structure that was seen at the outer edge of the SO envelope, with the 
velocity decreasing with the decreasing distance from the sour ce, is not seen h ere anymore. 
That PV structure is expected to be from the outflow walls (ILee et al.l 120001 ) . Thus, the 
ou ter edge of the SO envelope does not trace the inner edge of a pseudo-disk, as suggested 



m 



Paper t . As a result, the previous comparison of the inner part of the envelope to a 



collapsing magnetized rotating toroid of lAUen. Li. &: Shul ( 120031 ) and iMellon fc Lil (120081 ) is 
not appropriate anymore. 



5.2. Collapsing envelope 



The outer part of the envelope seems infalling towa rd the center with c o nstant specific 



angu l ar momentum, as found in other Class I sources (lOhashi et al.lll997al ; iMomose et al. 



19981 ; ILee &: Holl2005l ). The kinematics there is qualitatively consistent with that predicted 



in some infalling disklike enve lopes (or pseudo-disks) , either magnetized (jTerebev et al.lll984 



Krasnopolsky fc Konigll 120021 ) or nonmagnetized (see, e.g.. lNakamurall2000l ). where the ma- 
terial is infalling toward the center with its angular momentum. The infall velocity is ~ 80% 
of the rotation velocity at ~ 7000 AU at the edge of the envelope, and it drops to 40% of 
the rotation velocity at ~ 2000 AU at the transition radius. Therefore, in details, this part 
of the e nvelope is actually different from those collapsing envelop es found in other Class I 
sources (lOhashi et al.lll997ar IMomose et al.lll998t ILee fc Hoi 120051 ). where the infall velocity 
is larger than the rotation velocity. The HH 111 is probably in a later stage than those 
sources, with the motion already dominated by the rotation. The material in the envelope 
here is slowly spiraling inward and the centrifugal force progressively increases, so that the 
rotation can indeed become Keplerian at ~ 2000 AU. The infall rate in this part of the 
envelope is constant and given by 



M = 2.8m H / nv r r z dn w 4.2 x 10 -6 M© yr 



This infall rate is similar to that seen in other Class I sources (Ohashi et alll997a; Momose et al 



19981 ). Note that the infall rate here is a factor of 2 lower than that derived in lPaper t because 
of a factor of 2 lower in the infall velocity derived here. 
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An accretion shock is expected to occur at the transition radius, see e.g. jYorke fc Bodenheimer 

( 119991 ) . This accretio n shock has also been suggested to be present in the younger proto- 
stellar system L1157 (jVelusamy et al.l 120021 ) at a much smaller radius at ~ 400 AU at a 
much younger age in the Class phase. Here in HH 111, the isothermal sound speed in the 
envelope, which can be derived from the temperature, is ~ 0.15 km s" 1 at ~ 7000 AU and 
0.22 km s" 1 at 2000 AU. Thus, in the outer part of the envelope, the infall velocity is ~ 
1.1—1.4 times the isothermal sound speed, and is thus only slightly supersonic at that large 
distance. As a result, unlike that in LI 157, the accretion shock is expected to be weak at 
the transition radius and thus not easy to be detected. 



5.3. 2000 AU Sub-Keplerian disk? 

In the inner part of the envelope, where the rotation has become Keplerian, the infall 
velocity is uncertain and assumed to follow the same velocity profile as that in the outer part. 
It is ~ 40% of the rotation velocity. It increases from 1.4 to 2 times the isothermal sound 
speed from the transition radius to the vicinity of the source at ~ 240 AU (0"6). However, 
as mentioned earlier, the infall velocity could be overestimated and it could actually be 
subsonic as expected for a rotationally-supported disk with a Keplerian rotation. Therefore, 
the transition radius could indeed be the centrifugal radius, where the centrifugal force is 
balanced by the gravitational force. If that is the case, then the rotationally supported disk 
is quite large with a radius of ~ 2000 AU. A dis k with such large radius has been seen 



in the simulations of lYorke fc Bodenheimer! ( 119991 ). as a result of a gravitational collapse 



of a rotating spherical cloud. In their simulations, the angular momentum is transported 
outward via tidally induced gravitational torque. The HH 111 system can be compared to 
their model H, in which the cloud was set up to have a mass of 2 M & and a specific angular 
momentum of ~ 1333 AU km s _1 at the cloud edge. In that model, the cloud initially has 
a radius of 13333 AU and it collapses to form a quasi-static Keplerian disk surrounded by 
an infalling envelope. At the age of ~ 4 x 10 5 yr, the disk has grown up to ~ 1800 AU in 
radius with a central mass of 1.121 M & (see their Figures 2), similar to the HH 111 system. 
Further observations at higher resolution are needed to check if the infall velocity here can 
really be subsonic. 

Alternatively, this inner part of the envelope can be considered as a 2000 AU sub- 
Keplerian disk with a small infall motion. Assuming that both the infall and rotation 
motions are due to the gravitational force from the center, then the mass of the center 
source would be 1 + 0.5(v r /vk) 2 = 1.08 times that derived earlier, and the disk would be 
rotating at 1/VL08 ~ 96% of the Keplerian velocity. A similar sub-Keplerian disk with 
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a radius of ~ 2000 AU is also seen in the Class I source L1489 IRS, with a similar mass 
(~ 1.3 M Q ) of th e central source an d a similar ratio (~ 0.4) of the infall velocity to the 
rotation velocity (jBrinch et al.l 120071 ). A smaller rotating disk with a radius of ~ 600 AU 
is also seen in the Class I source IRAS 04365+2535, with a smaller mass of ~ 0.5 M Q for 



the central source (lOhashi et al.lll997bl ). This disk seems to have some small infall motion 
too Johashi et aD Il997bh and thus can also be considered as a sub-Keplerian disk. Note 
that the sizes of these sub-Keplerian disks seem to scale roughly linearly with the mass of 
the central source. These sub-Keplerian disks are expected to become fully rotationally- 
supported Keplerian disks closer to the source, where the infall velocity is expected to drop 
below the sound speed and thus become much smaller than the rotation velocity. Here in HH 
111, the rotationally-supported disk must have already formed around the source to launch 
the jet. 

The accretion rate in the disk can be assumed to be the same as the infall rate in 
the outer part of the envelope. From CO bullets along the j et axis, the mass-loss rate was 
estimated to be ~ 4xl0 -7 M & yr _1 in one side of the jet (ICernicharo fc Reipurthl Il996l ) 



and thus M~ ~ 8x10 7 M yr 1 in two sides. Therefore, the mass-loss r ate to accretion 



rate ratio is ~ 0.19 , comparable to that predicted in current MHD models (IShu et a. 



Pudritz et all l2007h . With a velocity of ~ 300-500 km s" 1 JCernicharo fc Reipurthl Il996h . 



2000; 



the jet has a mechanical luminosity of ^Mj vj ^ 6-16 L^, a bout 30-80% of the bolometric 
luminosity of the source, which is ~ 20 L & (lReipurthlll989l . corrected for the new distance 
of 400 pc). Assuming a stellar mass of ~ 1.3 M & and a stellar radius of ~ 3i? , the 
accretion luminosity is L acc = GM*M acc /R* ~ 58 L & . This luminosity is comparable to 
the sum of the mechanical luminosity of the jet and bolometric luminosity of the source, as 
expected if both luminosities are derived from gravitational potential energy released in the 
accretion. 



5.4. SO shocks in the disk? 



The SO emission is seen in the equatorial plane, with a mean radius of ~ 400 AU (1"). 
It is seen with a rotation similar to that seen in C 18 0. Since SO emission has been found to 



trace warm and shocked regions ( jCodella et al.l 120051 ; iLee et al.ll2006l . |2010| ). could the SO 
emission here trace a shock in and around the disk? The peak intensity of the SO emission 
is ~ 2 Jy beam -1 km s _1 in a beam of l!'3xl"l. Assuming that the emission i s optically 
thin and is from a region in LTE at a temperature of 100 K ( jCodella et al.l 120051 ). then the 
column density of SO is ~ 8 x 10 14 cm -2 . Since the emission is not well resolved, this value 
should be considered as a lower limit. Dividing this column density by that of H2 in the 
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model at 1", which is ~ 6.2xl0 22 cm -2 , the SO abundance is found to be ~ 1.3 x 10~ 8 , 
higher than that found in the ambient material, which is ~ l- 4x!0~ 9 (ICodella et al.ll2005l ). 
but lower than that found in the jet, which is ~ 2 xl0~ 6 (ILee et al.l |2010| ). Therefore, 
the SO emission here could trace a weak shock in the inner part of t he disk, such as the 
inner accretion shock on the disk surface (jYorke fc Bodenheimerl 119991 ). However, it is also 
possible that the SO emission traces the accretion shock between the sub-Keplerian disk and 
the rotationally-supported disk. 



5.5. Ejection? 

It has been proposed that the VLA 1 and 2 sources were first formed at the d ark ridge 



position and then ejected to the east and west, respectively (IReipurth et al.lll999l ). In our 
observations, the peaks of the C 18 emission shift from the dark ridge position to the east 
to the VLA 1 position from the low to high velocity. This shift of the emission peaks seems 
to support that the VLA 1 source was indeed formed at the dark ridge position and then 
eject ed to the east, as pro posed. The ejection was estimated to take place only ~ 3000 yrs 
ago ( Reipurth et al.|[l999l ). Therefore, the outer part of the envelope that is seen at low 



velocity does not have enough time to react to this ejection. It is aligned with the dark 
ridge likely because it is the one that produces the dark ridge with dust extinction. No dust 
emission is detected toward the dark ridge position because of the low density in the outer 
part of the envelope. On the other hand, the inner part of the envelope that is seen at high 
velocity and the disk have reacted to this ejection, carried away by the VLA 1 source, as 
seen in the observations. 



6. Conclusion 

The envelope and disks in the HH 111 prot ostellar system can be better studied with 



more complete uv coverage than that in iPaper t . The 1.33 mm continuum emission is now 



spatially resolved into two components, a bright one associated with the VLA 1 source and 
a much fainter one associated with the VLA 2 source, likely tracing the dusty putative disks 
around the two sources. The emission around the VLA 1 source is elongated in the equatorial 
plane with a radius of ~ 300 AU. The envelope is well seen in C 18 0, extending from ~ 60 
to 7000 AU out from the VLA 1 source, with the innermost part overlapping with the dusty 
disk. It has a differential rotation, with the outer part (~ 2000-7000 AU) better described by 
a rotation that has constant specific angular momentum and the inner part (~ 60-2000 AU) 
by a Keplerian rotation. The envelope seems to have some infall motion too. The amount 
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of infall motion, however, could be overestimated due to outflow contamination. The infall 
velocity to rotation velocity ratio is found to decrease from ~ 0.8 at the edge of the envelope 
to 0.4 at the transition radius and in the inner part of the envelope. Thus, the material in the 
outer part of the envelope seems to be slowly spiraling inward with its angular momentum 
and the rotation can indeed become Keplerian in the inner part. A compact SO emission is 
seen around the VLA 1 source with a radius of ~ 400 AU and it may trace a shock such as 
an (inner) accretion shock around the disk. It seems that the VLA 1 source was first formed 
in the dark ridge and then ejected to the east to its current position, as previously proposed. 

I thank the SMA staff for their efforts in running and maintaining the array. I also 
thank Sheng-Yuan Liu, Mike Cai, Ruben Krasnopolsky, and Yao-Yuan Mao for fruitful 
conversations. I also thank the referee for the valuable comments. 
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Fig. 1. — The 1.33 mm continuum emissi on contours overplotted on the HST NICMOS 
image adopted from iReipurth et al.l ( 119991 ). The gray and green contours are from the 
restored map and the CLEAN component map, respectively, of the continuum emission. 
The asterisk, cross, and "x" mark the positions of the VLA 1 and 2 sources, and the center 
of the dark ridge, respectively. The gray contour levels are 5a (1 — r n )/(l — r), where r = 1.3, 
n=l,2,3.., and a = 1.6 mJy beam -1 . The green contour levels start at 4%, with a step of 3%, 
of the peak value. The synthesized beam has a size of l'.'Axl'.'l. The blue and red arrows 
indicate the orientations of the blueshifted and redshifted parts of the jet, respectively. 
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Fig. 2. — C 18 and SO spectra toward the VLA 1 source averaged over an elliptical region of 
2"xl" in size with the major axis in the equatorial plane. The SO spectrum can be roughly 
described by a Gaussian line profile (solid curve), and its line center can be used here to 
define the systemic velocity. 
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Fig. 3. — Contours of blueshifted (blue) and redshifted (red) emission of C 18 and SO on the HST 
NICMOS image. The asterisk, cross, "x" , and the blue and red arrows have the same meanings as 
those in Figured! The dashed line indicates the equatorial plane perpendicular to the jet axis. The 
beam has a size of l"6xl"5 in (a) and l{'3xl{'l for the rest. In (a), a large beam is used to show 
the large-scale structure of the envelope. (Top two panels) Low-velocity blueshifted and redshifted 
emission of C 18 0. In (a), the contour levels start at 4a with a step of 4<r, where a = 0.032 Jy 
beam -1 km s _1 . In (b), the contour levels start at 3a with a step of 3a, where a = 0.03 Jy beam -1 
km s _1 . (Middle three panels) Low-, medium-, and high-velocity (i.e, LV, MV, and HV) blueshifted 
and redshifted emission of C 18 in the central region. The contour levels are the same as those 
in (b). (Bottom three panels) Same as those in the middle three panels, but for SO. The green 
contours in (e) and (/) are the same as those in Figure (H showing the CLEAN component map of 
the continuum emission. 
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Fig. 4. — Position-velocity (PV) diagrams in C 18 (black) and SO [gray in (d) and (e)] centered 
at the VLA 1 source. The gray boxes in the lower-right corners show the velocity and spatial 
resolutions of the PV diagrams, (a-d) PV diagrams cut perpendicular to the jet axis with increasing 
angular resolution toward the center. The resolutions are l"6xi''5 in (a), l"3xl"l in (6), and 
l''2xl{'0 in (c) and (d). The asterisks mark the data points that are used in Figure[5j Solid curves 
are derived from the rotation that has constant specific angular momentum. Dashed curves are 
derived from the Keplerian rotation. Here d = R/Rq, with R being the radial distance from the 
source and Ro =1" (400 AU). (e) PV diagrams cut along the jet axis. The resolution is the same as 
that in (6). The label "To VLA2" means that the SO emission there extends to the VLA 2 source. 
For C 18 0, the contours start at 2a and have a step of 2a, where a = 0.06 Jy beam -1 in (a), 0.045 
Jy beam -1 in (b) and (e), and 0.037 Jy beam -1 in (c) and (d). For SO, the contours start at 3a 
and have a step of 2a, where a = 0.045 Jy beam -1 in (d) and (e). 
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Fig. 5. — Rotation velocity and implied specific angular momentum plotted as a function of 
the radial distance from the source, with both corrected for the inclination angle. Open and 
filled circles are from the blueshifted and redshifted emission, respectively. Here d = R/R , 
with R being the radial distance from the source and R = 1" (400 AU). (c) and (d) are 
derived from (a) and (b), respectively, by shifting the center of the envelope by 0"15 to the 
south. 
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The best-fit parameters in our model for the C 18 envelope. 
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Fig. 7. — A simple model fitting to the integrated intensity map, spectrum, and PV diagrams 
in C 18 0. Black contours and spectrum are from the model. Gray contours with image and 
spectrum are from the observations, (a) shows the integrated intensity maps. The contours 
start at 3cr and have a step of 3a, where a = 0.06 Jy beam -1 km s -1 . (b) shows the spectra 
toward the VLA 1 source, (c) shows the PV diagrams cut perpendicular to the jet axis, as 
in Figure @J b). (d) shows the PV diagrams cut along the jet axis, as in Figure IIJe). 



